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Introduction 


It is well known that under some operating conditions, rocket engines (using solid or 
liquid fuels) exhibit unstable modes of operation that can lead to engine malfunction and 
shutdown. The sources of these instabilities are diverse and are dependent on fuel, 
chamber geometry and various upstream sources such as pumps, valves and injection 
mechanism. It is believed that combustion-acoustic instabilities occur when the acoustic 
energy increase due to the unsteady heat release of the flame is greater than the losses of 
acoustic energy from the system [1,2]. 

Giammar and Putnam [3] performed a comprehensive study of noise generated by gas- 
fired industrial burners and made several key observations; flow noise was sometimes 
more intense than combustion roar, which tended to have a characteristic frequency 
spectrum. Turbulence was amplified by the flame. The noise power varied directly with 
combustion intensity and also with the product of pressure drop and heat release rate. 
Karchmer [4] correlated the noise emitted from a turbofan jet engine with that in the 
combustion chamber. This is important, since it quantified how much of the noise from 
an engine originates in the combustor. 

A physical interpretation of the interchange of energy between sound waves and unsteady 
heat release rates was given by Rayleigh [5] for inviscid, linear perturbations. Bloxidge et 
al [6] extended Rayleigh’s criterion to describe the interaction of unsteady combustion 
with one-dimensional acoustic waves in a duct. Solutions to the mass, momentum and 
energy conservation equations in the pre- and post-flame zones were matched by making 
several assumptions about the combustion process. They concluded that changes in 
boundary conditions affect the energy balance of acoustic waves in the combustor. 
Abouseif et al [7] also solved the one-dimensional flow equations, but they used a one- 
step reaction to evaluate the unsteady heat release rate by relating it to temperature and 
velocity perturbations. Their analysis showed that oscillations arise from coupling 
between entropy waves produced at the flame and pressure waves originating from the 
nozzle. Yang and Culick [8] assumed a thin flame sheet, which is distorted by velocity 
and pressure oscillations. Conservation equations were expressed in integral fonn and 
solutions for the acoustic wave equations and complex frequencies were obtained. The 
imaginary part of the frequency indicated stability regions of the flame. 

Activation energy asymptotics together with a one-step reaction were used by McIntosh 
[9] to study the effects of acoustic forcing and feedback on unsteady, one-dimensional 
flames. He found that the flame stability was altered by the upstream acoustic feedback. 
Shyy et al [10] used a high-accuracy TVD scheme to simulate unsteady, one-dimensional 
longitudinal, combustion instabilities. However, numerical diffusion was not completely 
eliminated. Recently, Prasad [11] investigated numerically the interactions of pressure 
perturbations with premixed flames. He used complex chemistry to study responses of 
pressure perturbations in one-dimensional combustors. His results indicated that reflected 
and transmitted waves differed significantly from incident waves. 
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In the present paper, an extension of the work perfonned by Frendi [12] is carried-out. In 
particular, an attempt is made at understanding the onset of combustion instabilities as 
well as the use of passive devices to control them. Following the systematic experimental 
investigation of Laudien et al. [13], who studied the design of acoustic cavities, a 
numerical investigation is performed. The remainder of the paper is organized as follows; 
a brief description of the mathematical model is given in the next section followed by the 
numerical techniques used to solve it. Details of the results are then discussed followed 
by the concluding remarks. 

Mathematical Model 


The model used to describe this problem is based on the nonlinear Euler equations, which 
can be written as 


au fF, 

at ax i 
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and E = c v T + / 2 v l v i . The closure is obtained through the equation of state for an ideal 
gas, p = pRT . For a two-dimensional axi-symmetric model the definition of the vector B 
results from the transformation to cylindrical coordinates. If xi is the axis of symmetry, 
the vector B can be written as 
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Method of Solution 


The numerical methodology employed to solve Equation (1) can be described as a 
combination of a temporal discretization scheme, a spatial discretization scheme and, a 
discontinuity-capturing scheme. In the time discretization scheme the value of U is 
sought given its value at a previous time by means of a Taylor series expansion. The 
spatial discretization is done using the Galerkin finite element method in which the 
integrated weighted-residual is minimized. To resolve discontinuities in the flow field, a 
convective flux correction term is devised. A second order accuracy is achieved in time 
while a second order or higher accuracy is achieved in space. Details of the numerical 
scheme can be found in [14, 15]. The boundary conditions used are rigid slip wall along 
the chamber walls and characteristic boundaries at the inlet (upstream boundary) and exit 
(downstream boundary). When acoustic disturbances are introduced in the chamber, the 
inlet pressure is specified as 
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p = p 0 (\ + £sm(a>t)) or 


p = p 0 (\ + sR(t)) 


( 4 ) 


with s being the excitation amplitude, co = 2rf the frequency and R(t ) a random 
number. In Equation (4), p 0 is the chamber pressure. The other inlet quantities are 
specified using the method of characteristics. 

Results and Discussions 


The numerical investigation carried-out in this paper follows the detailed experimental 
investigation performed by Laudien et al. [13]. In particular, acoustic tests are performed 
in an axi-symmetric combustion chamber having the shape shown on Figure 1. A 
loudspeaker is placed at the center of the inlet boundary and emits plane harmonic or 
random acoustic waves. As described in section III, a finite element technique is used. 
Figure 1 also shows a typical grid used. A denser grid is used near the acoustic source. 
The conditions in the chamber are those of air at T = 300 K and P = I atm. 



Figure 1: Typical unstructured grid used. 

At first, the natural frequencies of the chamber are determined. To this end a random 
acoustic disturbance is introduced by the loudspeaker and the pressure fluctuations are 
monitored at several points in the chamber. Figure 2 shows the sound pressure level 
(SPF) as a function of the frequency at a point along the top chamber wall. Several peaks 
corresponding to the various chamber modes are obtained. In order to determine the 
mode shapes, one can excite the individual modes by introducing a plane, hannonic 


XV-4 


acoustic wave at the given frequency. Following the determination of the natural 
frequencies of the chamber, a frequency of 1880 Hz, corresponding to a chamber mode, 
is selected for extensive tests. The first such test is to obtain the chamber response to a 
high (i.e. nonlinear) and low (i.e. linear) excitation amplitude. Figures 3(a)-(b) show the 
time histories of the pressure at an observation point along the top wall. 



Figure 2: Frequency response of the combustion chamber to a random excitation. 

Figure 3(a) shows rapid growth of the pressure oscillations and a slowly varying mean 
pressure, whereas Figure 3(b) shows a slow growth of the oscillations around a zero mean 
pressure. The corresponding power spectra are shown on Figs. 4(a)-(b). For the high 
excitation amplitude case, s = 0.5 , strong hannonics are shown to exist along with the 
fundamental frequency. It is believed that the nonlinear interaction between the fundamental 
and its harmonics is responsible in part for the rapid growth of the fluctuations. This point is 
confirmed by Fig. 4(b) which shows weak harmonics and therefore a quasi-linear behavior 
and weaker growth as shown by Fig. 3(a). 
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Figure 4: Sound pressure level at a point along the top wall (a) s - 0.5 , (b) s = 0.001 . 

In order to assess the level of damping of a given mode by the chamber, the loudspeaker is 
turned on and then off and the pressure oscillations in the chamber are monitored at a given 
point. Figure 6 shows the damping level for two different modes 1880 Hz (Fig. 5(a)) and 512 
Hz (Fig. 5(b)). In addition to the pressure oscillations, a decay curve described by Xe~ n<u ' is 
also shown, where 77 is the decay factor. The Figure shows that the mode corresponding to 
1880 Hz decays faster than the mode corresponding to 512 Hz. 
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Figure 5: Decay rate of the pressure oscillations for two different chamber modes 
corresponding to (a) 1880 Hz (b) 512 Hz. 


Baffles 


Following the preliminary studies presented above, the effect of adding a baffle in the radial 
direction of the inlet boundary, Figure 6, on the pressure oscillations in the chamber, is 
investigated. Figure 7 shows the pressure oscillations at the observation point along the top 
wall obtained with and without the baffle. The Figure shows a significant reduction of the 
oscillation amplitude of the pressure. The excitation amplitude used is s — 0. 1 . 
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Figure 6: Sketch of the inlet boundary with a circular baffle; (a) front view, (b) side view 
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Figure 7: Effect of a radial baffle on the pressure oscillations in the chamber. 


Figure 8 shows the corresponding power spectra. This figure answers a critical question that 
can be asked based on Figure 7, which is: where did the excitation energy go? The answer 
based on Figure 8 is: the energy is redistributed on a broader frequency spectrum. The reason 
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for this redistribution is that the excitation frequency, 1880 Hz, is no-longer a fundamental 
mode of the chamber, therefore there is no resonance. 



Figure 8: Power spectra of the pressure oscillations in the chamber with and without a baffle. 


The effect of the baffle on the mode corresponding to the 512 Hz frequency is shown on 
Figs. 9(a)-(b). Figure 9(a) shows that after the excitation source is turned-off, Time > 30, a 
weak decay of the pressure oscillations is obtained without using the baffle. When a baffle is 
used, Fig. 9(b), the decay rate is larger as indicated by the steeper exponential curve. 
Therefore, one can conclude that even though the mode corresponding to the 512 Hz 
frequency is not suppressed by the baffle, its damping is increased. This results is in 
agreement with the experimental measurements of Laudien et al. [13]. The overall sound 
pressure level in the chamber based on the root-mean-square (rms) pressure has also dropped 
from 170 dB to 156.5 dB. 

Helmholtz Resonators 


Helmholtz resonators, known also as acoustic cavities, are very popular passive control 
devices. Their use is wide spread; from controlling instability waves to noise reduction 
devices such as liners. A Helmholtz resonator is composed of an inlet of area S and length l 
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Time 

Figure 9: Time history of the chamber pressure along the top wall with the excitation source 
turned-on and then off (Time > 30), (a) without a baffle, (b) with a baffle. 


followed by a cavity of volume V, as shown in Fig. 10. For such a resonator, the resonant 
frequency is given by 

fo = 
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V s 



where c is the speed of sound in the chamber and Al ~ 0.85 d with d being the inlet 
diameter. The Helmholtz resonator used is tuned to a resonant frequency of 1880 Hz, which 
corresponds to a chamber mode. When all the corners are sharp. Figure 1 1 shows that the 
presence of the resonator has reduced the pressure fluctuations in the chamber significantly. 
Similar to the baffle case, the power spectra show the presence of more frequencies when the 
resonator is used, supporting the argument that energy has been diverted from the 
fundamental to other frequencies. The overall sound pressure level in the chamber based on 
the rms-pressure has also dropped from 170 dB to 154 dB (2.5 dBs better than the baffle). 
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Figure 11: Effect of a Helmholtz resonator on the pressure oscillations in the combustion 

chamber. 
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Figure 12: Power spectra of the pressure oscillations in the chamber with and without a 

resonator. 


The effect of smoothing the various corners of the resonator on the overall sound pressure 
level in the chamber is shown on Figure 13. The figure shows that the best result is obtained 
when both inlet corners (referred to as throat corners on the Figure) are smoothed. Smoothing 
the other interior resonators comers did not achieve any additional improvement, as shown 
by the figure. Notice that “inlet corner” on the figure refers to the corner facing the 
combustion chamber. 
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Figure 13: Effect of smoothing the resonator comers on the overall sound pressure level. 


According to the experimental results of Laudien et al. [13], the resonator becomes less 
effective at higher noise levels in the chamber. This result is confirmed by our computations, 
which show less-drop in dB-level, Figure 14, even with the best configuration of Fig. 13. 
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Figure 14: Effect of a Helmholtz resonator on the overall sound pressure level of a 

combustion chamber at high noise levels. 

The problem in designing a Helmholtz resonator is the required knowledge of the tuning or 
resonant frequency. Therefore, one question we needed to answer is how effective a 
resonator is when the design frequency is slightly off? This is of practical importance since in 
an actual combustion experiment the frequency to be controlled is not known accurately. To 
answer this question, a resonator tuned to a frequency of 1800 Hz was used to control a mode 
of frequency 1880 Hz. Figure 15 shows the overall sound pressure level in the chamber with 
an “On-Design” and an “Off-Design” resonator. The “Off-Design” resonator achieves 13.5 
dB less reduction than the “On-Design” one. This is significant since the “Off-Design” 
frequency is only 80 Hz off the design value. 


XV- 15 




Figure 15: Effect of resonator design frequency on the overall sound pressure level in the 

combustion chamber. 


)uarter-Waves 


Another passive control device commonly used is a quarter-wave. It is somewhat similar to a 
Helmholtz resonator except for the absence of the cavity part, Figure 16. The resonant 
frequency of a quarter-wave is given by 


4 (L + M) 


with c being the speed of sound and M ~ 0.85 d with d the quarter-wave diameter. 
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Figure 16: Typical geometry of a quarter-wave. 


In practice, several shapes of quarter-waves are used. In this study two of these shapes will 
be investigated; an L-shaped quarter-wave and a straight quarter-wave. In addition, the effect 
of smoothing the comers will be analyzed. Figure 17 shows the time history of the chamber 
pressure with and without a tuned quarter-wave. A significant reduction in pressure 
oscillations is achieved when using a quarter- wave. Figure 18 shows the corresponding 
power spectra. Similar to the baffle and resonator, the quarter-wave effect is to detune the 
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Figure 17: Effect of a quarter-wave on the pressure oscillations in the combustion chamber. 
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Figure 18: Power spectra of the pressure oscillations in the chamber with and without a 

quarter- wave. 



Figure 19: Effect of L-shaped and Straight quarter-waves with different corner geometries 

on the overall sound pressure level in the chamber. 

chamber and spread the energy to other frequencies. A comparison of an L-shaped quarter- 
wave to a straight one with various comer geometries is shown on Figure 19. For all corner 
geometries, the straight quarter-wave achieves more reduction in SPL than the L-shaped one. 
The figure also shows that smoothing the inlet comers has an adverse effect on the overall 
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SPL in the chamber for both quarter- wave geometries. For the straight quarter- wave, a 
narrow inlet resulted in a higher SPL, as expected. 

Figure 20 shows a comparison of overall SPL in the combustion chamber for all the passive 
devices tested in this study. The figure shows that the Helmholtz resonator with a smooth 
inlet gives the best results. 



Figure 20: Comparison of overall sound pressure level in the chamber obtained with the 

various passive devices. 

Figures 21-24 show the instantaneous pressure contours in the chamber obtained with the 
various passive devices. A different pattern can be observed for each device. 
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Figure 21: Instantaneous pressure contours in the combustion chamber obtained without a 

passive device. 



Figure 22: Instantaneous pressure contours in the combustion chamber obtained with a 

baffle. 
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Figure 23: Instantaneous pressure contours in the combustion chamber obtained with a 

Helmholtz resonator. 



Figure 24: Instantaneous pressure contours in the combustion chamber obtained with a 

quarter- wave. 
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Concluding Remarks 


Numerical experiments have been carried-out to study the effectiveness of various 
passive devices to control combustion instabilities. The results show that, for the given 
mode studied, a Helmholtz resonator with a smooth inlet achieves the best control results. 
However, at high sound pressure levels in the chamber, the resonator becomes less 
effective. For the case of a quarter-wave, it is found that a straight shape performs better 
than an L-shaped quarter- wave. In addition, smoothing the comers of a quarter- wave inlet 
may have an adverse effect on the overall sound pressure level in the chamber. The 
computational results presented here are in good qualitative agreement with the 
experimental measurements of Laudien et al. [13]. 
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